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Resumo 

A utilização de imagens de alta resolução espacial para fins de planejamento urbano, na determinação de parâmetros relacionados ao uso do 
solo como a fração de áreas impermeáveis correlacionando a densidade habitacional, é atualmente uma ferramenta presente e indispensável 
aos tomadores de decisão. Uma pequena revisão de literatura aborda a estreita relação entre a densidade populacional e a impermeabilidade 
do solo dada a ocupação urbana. Essa estreita relação invoca ao planejamento urbano elementos reguladores da expansão urbana, já que a 
intensa impermeabilização do solo afeta diretamente o ciclo hidrológico pela redução da infiltração das águas pluviais e consequente aumento 
do escoamento superficial, proporcionando, nos períodos chuvosos, alagamentos frequentes e inundações. Deste modo, a quantificação de 
áreas impermeáveis permite aos gestores públicos definir ações relacionadas à infraestrutura como transporte, saneamento ambiental, saúde 
e educação. Utilizando imagem de alta resolução espacial (lm) do satélite QuickBird II, obtiveram-se para 12 bairros da cidade de Porto 
Alegre-RS a porcentagem de áreas impermeáveis e sua relação com a densidade habitacional por meio de dados do censo demográfico. Para 
verificação da classificação rígida utilizada, o estudo apoiou-se na obtenção do acerto global e do coeficiente de kappa, obtendo-se valores 
respectivamente de 70% e 0,57, considerados como satisfatórios para o tipo de classificação utilizada. Pelas informações obtidas, percebe-se 
que apesar do decrescimento populacional ocorrente nas regiões centrais e o avanço sobre a periferia da cidade, não houve impedimento 
para a densificação no núcleo central com aumento de áreas impermeáveis. 

Palavras-chave: : geoprocessamento, densidade habitacional, classificação rígida 


Abstract 

The use of high spatial resolution satellite image for urban planning to determination of parameters related to land use as the fraction of 
waterproof areas correlating the housing density is currently a present and indispensable tool for decision makers. A short review of literature 
discusses the close relationship between population density and impermeability of the soil due to urban occupation. This close relationship 
refers to urban planning, regulatory elements of urban sprawl, as the intense waterproof surfaces directly affects the hydrological cycle by 
reducing the infiltration of rainwater and consequent increase in runoff, providing frequent flooding in the rainy season and floods. Thus, 
the quantification of waterproof areas, allows managers to actions related to public infrastructure such as transport, sanitation, health and 
education. Using high spatial resolution image (lm) of QuickBird II was obtained for 12 neighborhoods of the city of Porto Alegre-RS, the 
percentage of waterproof areas and its relationship with housing density by census data. To verify the used hard classification, the study was 
supported in obtaining the overall arrangement and the kappa factor, obtaining values of respectively 70 and 0.57% considered as satisfactory 
for the classification type used. From the information obtained it is clear that despite the population decline occurring in the central regions 
and advance on the outskirts of the city, there was no impediment to the densification of the central nucleus with increased waterproof areas 

Keywords: Keywords: geoprocessing, housing density, hard classification. 


Recebido em: 18.05.12, Revisado em 21.02.13, Aceito em 29.03.13 


34 


Menezes Filho | Estimativa das áreas.... 


I. INTRODUÇÃO 

A utilização de imagens de alta resolução 
para determinação de áreas impermeáveis em áreas 
urbanas torna-se uma prática comum nos seto- 
res de planejamento urbano. Vantagens quanto à 
precisão dos resultados e tempo envolvido para a 
classificação de imagens na determinação de áreas 
impermeáveis baseadas em imagem de satélite são 
relatadas desde a década de 70, conforme trabalho 
de Ragan e Jackson (1975). 

A necessidade da obtenção da fração de áreas 
impermeáveis se relaciona à densidade habitacional 
e, assim, constitui-se ponto de apoio ao ordena- 
mento territorial e previsão de cenários futuros de 
crescimento urbano. Do mesmo modo, a parcela 
de áreas impermeáveis é utilizada como parâmetro 
de entrada de modelos hidrológicos utilizados para 
previsão de cheias urbanas. 

Consoante Jensen (2009), a inferência popu- 
lacional pode ser realizada a nível local, regional e 
nacional, tendo como base: a contagem individual 
de unidades domiciliares - sendo a mais precisa, 
contudo inviável para escalas regional e/ou nacional 


devido ao tempo e custo operacional para realização; 
a medição da extensão de áreas urbanizadas e por 
meio de estimativas derivadas da estreita relação 
da área urbana com a densidade populacional que 
é foco deste trabalho. 

Dentre os trabalhos pioneiros a nível inter- 
nacional sobre o assunto, destaca-se Stankowski 
(1972) pela proposição do equacionamento baseado 
em levantamentos de campo para 21 municipali- 
dades de New Jersey (EUA). No mesmo trabalho, 
demonstraram-se os diversos tipos de curvas base- 
adas nas características do uso do solo tais como: 
residencial, comercial e industrial, apresentando, as 
modificações de áreas urbanas e suburbanas pelo 
adensamento populacional e sua importância para 
o planejamento e estudos hidrológicos. 

Gluck e McCuen (1975), visando desen- 
volver equações pela rapidez na obtenção de áreas 
impermeáveis em contraposição ao tempo gasto 
na análise de fotografias aéreas, acrescentaram 
também outros parâmetros como: densidade ocu- 
pacional, densidade de moradias e a distância ao 
centro comercial a fim de refinar o estudo referente 
a Washington, E.U.A. 


Impermeável 



Figura 1. Modelo V-I-S (Vegetação-Impermeabilidade-Solo). Adaptado (RIDD,1995). 
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Huber e Dickinson (1992) citam Heaney et. 
al com uma coletânea das equações de regressão, 
obtidas de estudos de outros autores. Os mesmos 
autores definem a equação proposta por Stankovisky 
(1972) como a mais representativa pela dimensão 
do estudo realizado com 567 municipalidades e 
utilizada no modelo de gerenciamento de águas 
pluviais SWMM ( Stormwater Management Modet) 
da agência ambiental norte-americana ( EPA -US 
Environmental Protection Agency). 

Além de equações propostas, outras meto- 
dologias foram elaboradas para correlacionar área 
impermeável e densidade habitacional como o 
trabalho de Ridd (1995) que propôs um modelo 
denominado V-I-S (Vegetação - Impermeável - Solo) 
(Figura 1) por meio de projeto piloto com 770 pontos 
da cidade de Salt Lake nos E.U.A, representando 
a composição de cobertura urbana. zz 

Relacionado aos impactos urbanos no ciclo 
hidrológico, o modelo V-I-S demonstra claramente 
as relações do aumento de áreas impermeáveis com 
consequente redução da vegetação e decréscimo da 
porcentagem de solo exposto. Um dos primeiros 
trabalhos a nível nacional foi o realizado por Motta 
e Tucci (1984), utilizando ortofotos para bacias de 
Porto Alegre, no intuito de determinar a relação 
área impermeável versus densidade habitacional. 
Posteriormente, Tucci et. al. (1989) desenvolve- 
ram o mesmo estudo para onze bacias urbanas na 
região Metropolitana de São Paulo. Tais estudos 
foram conduzidos a grandes bacias urbanas porém 
restritos à aplicação das equações para condições de 
relevos distintas daquelas estudadas, como regiões 
de alta declividade e morros. Tal estudo contemplou 
ainda a estimativa para as cidades de São Paulo e 
Curitiba obtendo uma tendência geral da relação 
população-impermeabilidade (Figura 2). 


O surgimento de imagens comerciais de 
alta resolução, como a do satélite QuickBird II no 
início dos anos 2000, subsidiou para áreas urba- 
nas a identificação de alvos com maior precisão 
sobre as imagens anteriormente disponíveis como 
a LANDSAT ( LandRemote Sensing Satellitè), Cbers 
(China-Brazil Earth-Resources Satellitè), Aster ( Advan- 
ced Spaceborn Thermal Emission ) e Spot ( Système 
Pour VObservation dela Terrè). O tipo de classificação 
rígida baseada na Máxima Verossimilhança tam- 
bém se tornou frequente, conduzindo a resultados 
às vezes satisfatórios, outras vezes a resultados 
irreais, devido à confusão espectral caracterizada 
pelo a pixel mistura" como observado em diversos 
trabalhos (BERNARDI et al., 2007; CENTENO 
et al., 2003; PINHO e KUX; 2004). 

Campana e Tucci (1994) puderam utilizar 
imagens LANDSAT para determinação de áreas 
impermeáveis para a cidade de Porto Alegre e 
contemplaram uma área de 140 km 2 adotando uma 
abordagem fuzzy. A mesma imagem foi utilizada por 
Alves (2004) para determinação da mesma relação 
bem como da evolução de áreas impermeáveis na 
Bacia do Arroio Dilúvio na cidade de Porto Alegre. 

As limitações da relação área impermeável 
versus densidade habitacional compreendem: (a) 
áreas predominantemente comerciais e industriais 
não mostram a mesma tendência; (b) a topografia 
pode alterar a ocupação devido à dificuldade de 
ocupação (grande declividade do terreno); (c) áreas 
pequenas podem estar distorcidas por misto de 
edifícios com áreas verdes. 

Este trabalho buscou abordar a estimativa 
das áreas impermeáveis, baseada na classificação 
rígida de imagem de alta resolução espacial a fim 
de demonstrar sua estreita relação com a densidade 
habitacional e o planejamento urbano. 



Figura 2. Tendência geral da relação ocupação-impermeabilidade (CAMPANA e TUCCI, 1994). 


36 


Menezes Filho | Estimativa das áreas.... 


2. Materiais e métodos 

2. 1 . Classificação de imagens digitais 

A classificação de imagens multiespectrais 
consiste na identificação de diferentes alvos na 
superfí- cie terrestre através do comportamento 
espectral dos mesmos. Este processo pode par- 
tir do próprio usuário na identificação visual e 
no agrupamento de tais alvos em classes, sendo 
denominada de classificação supervisionada. O pro- 
cesso de classificação pode ser não-supervisionada 
ficando a cargo do software utilizado agrupar tais 
alvos a semelhantes comportamentos espectrais 
( “clustering”). 

A classificação supervisionada pode ainda 
ser subdividida em rígida ou relativa. Quanto à 
primeira classificação, destaca-se o método da 
máxima veros- similhança, um dos mais utilizados. 
Este método se baseia na teoria de Bayes. Conhe- 
cendo as informações do conjunto de amostras de 
treinamento de cada classe, este comando associa 
através da média, da variância e covariância das 
assinaturas espectrais, a possibilidade de um pixel 
pertencer a uma determinada classe (CLARK E 
LABS, 2006). Esse processo não contempla para 
o pixel a pertinência em diversas classes e atribui 
a este, como o nome diz, unicamente a classe de 
maior preponderância, independentemente se para 
o mesmo pixel houverem outras classes existentes 
em menor porcentagem. 

Na classificação relativa, baseada na Teoria 
dos Conjuntos Fuzzy, proposta por Zadeh (1965), 
interpreta-se como uma questão de grau a pertinên- 
cia de um objeto a uma classe. Nesta classificação, 
são geradas imagens que representam a probabili- 
dade dos pixels existentes na imagem pertencerem 
a cada classe, num intervalo de 0 a 1 . 

2.2 Area de estudo 

A área de estudo compreendeu 12 bairros da 
cidade de Porto Alegre-RS. Estes bairros possuem 
características amostrais distintas com relação à 
topografia e ocupação descritas na tabela 1 . 

Os dados demográficos foram obtidos no 
SIDRA, um banco de dados agregados que dispo- 
nibiliza conforme orientação do usuário na página 
principal do IBGE - Instituto Brasileiro de Geo- 
gráfica e Estatística, informações sobre contagem 
demográfica, por estados, cidades, municípios e 
bairros brasileiros consoante o censo demográ- 
fico do ano de 2000. Paralelamente, obtiveram-se 
informações sobre a área dos bairros na página da 


internet Observa POA (http://www2.portoalegre. 
rs.gov.br/observatorio/). 

Tabela 1 . Características dos bairros de estudo 


Bairro 

DH (hab/ha) 

Ocupação 

Bom Fim 

298,71 

Residencial 

Cidade Baixa 

210,56 

Residencial 

Centro 

153,10 

Residencial 

Santana 

142,42 

Residencial 

Menino Deus 

137,57 

Residencial 

Floresta 

89,47 

Com.Ind.Serv. 

Hípica 

23,18 

Residencial 

Guarujá 

20,55 

Residencial 

Serraria 

16,79 

Residencial 

Chapéu do Sol 

5.24 

Residencial 

Ponta Grossa 

3,53 

Residencial 

Lajeado 

1,26 

Residencial 


Fonte: Secretaria Municipal de Planejamento, SMP 


Buscou-se selecionar bairros representativos 
com características demográficas distintas para 
o estudo. Na Figura 3, encontram-se os bairros 
selecionados. 

De posse destes dados, partiu-se para a uti- 
lização de geoprocessamento por meio do software 
IDRISI Andes. A imagem utilizada para realização 
do trabalho foi a do satélite QuickBird //, imagem 
fusionada (pancromática mais as bandas do visível) 
com resolução espacial de lm, obtidas no período 
de março de 2002 a março de 2003. 

Com os bairros escolhidos, gerou-se então 
para cada bairro uma máscara para classificação 
individual (Figuras 4 e 5). Prosseguindo o trabalho e 
como etapa para classificação, foi necessário , para 
cada bairro, gerar as bandas que formam a imagem 
fusionada, utilizando o comando SEPARAT e, assim, 
gerando para cada bairro três arquivos compostos 
pelas bandas Vermelho’, Verde’ e ‘azul’ (RGB). 

Para aceitação da classificação, utilizou-se 
uma imagem de controle (100m x 100m) (Figura 6) 
referente a uma área pertencente a um dos bairros 
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Densidade Habitacional em bairros de Porto Alegre 


r/3 
r i 


Lageado 
Ponta Grossa 
Chapéu do Sol 

Serraria 

Guatvjà 

Hípica 

. B Floresta 

cc 

PQ Menino Deus 

Santana 

Centro 

Cidade Baixa 

Bom Fim 



100 


200 


300 


400 



Figura 3. Densidade Habitacional dos bairros selecionados 



Figura 4. Máscara gerada para o Bairro Lageado 



Figura 5. Máscara gerada para o Bairro Centro 
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de estudo com características de áreas impermeáveis 
e permeáveis comuns aos outros bairros de estudo. 

Utilizou-se como índice da qualidade da 
classificação o percentual global de acerto obtido por 
meio da matriz de erro ou confusão e o índice lcappa. 
Este índice busca representar o grau de precisão 
temática entre a imagem classificada e a imagem 
referência tomada como verdade, variando de 0 a 
1 . Quanto mais próximo 1 , têm-se resultados ditos 
como excelentes e próximos a 0 de baixa precisão 


Tabela 2 - Alvos selecionados para classificação 


Classes 

/ 

Area 

permeável 

/ 

Area 

Impermeável 

r 

Agua 


X 

Asfalto 


X 

Cobertura cerâmica 


X 

Cobertura concreto 


X 

Cobertura metálica 


X 

Vegetação 

X 


Solo exposto 

X 



A imagem de controle foi classificada visua- 
lmente, definindo-se duas classes principais, uma 


Wihilowfroiii Stintdiiti coinposite c: 1618 r: 1567 to c: 1717 r: 1666 



Figura 6. Imagem de controle para classificação visual 

permeável, englobando sub-classes como vegetação 
e solo exposto e outra classe impermeável, composta 
pelas sub-classes: asfalto, cobertura cerâmica, cober- 
tura concreto, cobertura metálica e água (Tabela 2) 
Realizou-se o cruzamento da imagem de controle 
classificada visualmente, admitindo-se esta como 
verdadeira, com a gerada pela classificação auto- 
mática, obtendo-se então a matriz de confusão e o 
cálculo do índice Kappa. Conforme Alvez (2004), 
este cruzamento permite avaliar o grau de precisão 
do que se imagina na realidade e a imagem classi- 
ficada automaticamente pelo software. 

Além disso, com a matriz de confusão gerada, 
pode-se ainda optar pelo maior refinamento ou 
escolha de um maior número de amostras e até 
mesmo a exclusão de amostras que geraram maior 
confusão espectral. O processo encontra-se ilustrado 
na Figura 7. 


IMAGEM DE CONTROLE 




1 




llfewlm 


«gfi 




Composição RGB 
BAIRRO 


SEPARATE 

(RGB) 



3and 1 extracted from corta santana 


3and 2 extracted from corte santana 



3and 3 extracted from corte santana 


CLASSIFICAÇÃO VISUAL 




Maximum Llkellhood Classification 


-> 



CROSSTAB 


Croi*-CU»»tfie«ion : tremamemovúual | cone maxllke 




Figura 7. Esquema ilustrativo do processo para avaliação da classificação automática 
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Figura 8. Classificação por Máxima Verossimilhança 




Após essa etapa, realizou-se a classificação 
automática para os bairros em estudo por Máxima 
Verossimilhança, utilizando as classes já citadas, 
posteriormente reclassificadas como permeáveis 
e impermeáveis (Figura 8), com os valores 1 e 2 
respectivamente. 

3. Resultados e discussão 

A classificação automática, utilizando 
Máxima Verossimilhança disposto no software 
Idrisi Andes , apresentou a mesma característica de 
trabalhos já desenvolvidos quanto à existência de 
confusão entre classes, como entre solo exposto e 
cobertura cerâmica e asfalto e vegetação. 

Os maiores erros de inclusão e omissão 
ocorreram para as classes asfalto, concreto e solo 
exposto, ocasionados pela confusão espectral e pela 
dificuldade apesar da alta resolução espacial, da 
definição das classes de treinamento pela presença 
de “pixel mistura”. Determinou-se, então, apriori, o 
erro total da classificação empregada na imagem de 
controle, obtendo-se o valor de 42%, uma exatidão 
de 58% e um índice kappa de 0,47, significando 
este resultado uma precisão média. Para contor- 
nar o problema intrínseco à classificação rígida de 
imagens de alta resolução espacial, aumentou-se 
a quantidade de amostras de treinamento, selecio- 
nando amostras representativas de solo exposto e 
asfalto, classes com maiores erros de omissão. 

Além da classe água representar 2% dos 
pixels e não estar presente na maioria dos bairros, 
com erro de inclusão de 99%, tal classe foi desconsi- 
derada. Do mesmo modo, constatou-se na imagem 
utilizada a grande presença de sombras de prédios 
sobre prédios e sobre o asfalto, claramente áreas 


impermeáveis nos bairros: Cidade Baixa, Menino 
Deus e Floresta. Ao não considerar a água como 
classe para os bairros onde não se detectou detectou- 
-se a mesma bem como a retirada de sombra da 
imagem de controle e a maior amostragem para 
as classes de maior confusão espectral como solo 
exposto e asfalto, a exatidão da classificação foi de 
aproximadamente 70% um índice kappa obtido de 
0,57, considerando-se satisfatória a qualidade da 
classificação. 


r 

Tabela 3. Valores de Area Impermeável nos bairros 


Bairro 

HD (hab/ha) 

Al (%) 

Bom Fim 

298,71 

71,58 

Cidade Baixa 

210,56 

73,12 

Centro 

153,10 

78,66 

Santana 

142,42 

78,33 

Menino Deus 

137,57 

69,70 

Floresta 

89,47 

68,04 

Hípica 

23,18 

20,82 

Guarujá 

20,55 

40,88 

Serraria 

16,79 

22,29 

Chapéu do Sol 

5.24 

15,60 

Ponta Grossa 

3,53 

13,25 

Lagrado 

1,26 

18,15 
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Destaca-se que a retirada das classes água 
e sombra se aplicou somente nos bairros onde tais 
conjuntos não foram relevantes em relação às carac- 
terísticas da área do bairro em estudo. Reafirma-se 
que, por tratar de área urbana, tais classes deverão 
compor a classificação por constituírem elementos 
presentes dentro dos alvos urbanos. 

Partiu-se, então, para a classificação auto- 
mática dos bairros em estudo considerando as 
classes já citadas. Determinou-se, então, a área 
impermeável para cada bairro. Os valores para as 
áreas impermeáveis são apresentados na Tabela 3. 

O bairro Floresta, por ter a ocupação refe- 
rente a serviços, comercial e industrial, apresentou-se 
de maneira clara pela baixa densidade e alta fração 
de áreas impermeáveis proporcionalmente ao con- 
junto de bairros em estudo (89,47 hab/ha, 68,04%). 

Analisando a Figura 9, fica evidente, con- 
soante Tucci (2007) que a relação cresce até 120 
hab/ha quando a partir desde momento o efeito 
da verticalização reduz a área impermeável pro- 


porcional próximo de 70% da área total. 

Segundo Alves (2004), a bacia do dilúvio 
nas quais os bairros foram amostrados apresentou 
um crescimento de 1987 a 2000 de 10% na taxa 
de urbanização, evidenciando uma alteração das 
características populacionais e de áreas imperme- 
áveis na área de estudo. 

Entretanto, analisando a evolução da popu- 
lação para os bairros em estudo, houve um decres- 
cimento populacional para bairros mais antigos e 
com ocupação consolidada como Bom Fim, Cidade 
Baixa, Centro e Santana. Tal análise demonstrou 
que de certa forma houve uma densificação das 
áreas existentes e possivelmente o aumento de 
áreas impermeáveis, justificando o incremento de 
10% na taxa de urbanização. 

O único crescimento verificado foi para 
bairros periféricos distantes do núcleo central den- 
tre eles: Guarujá, Serraria e Hípica com ocupação 
residencial/rural como observa-se na Tabela 4. 


Densidade Habitacional x Área Impermeável 



Densidade Habitacional (hab/ha) 


Figura 9. Relação densidade habitacional x área impermeável para bairros da cidade de Porto Alegre - RS. 
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Tabela 4. Variação populacional entre censos 


Bairro 

Pop. 1991 

Pop. 

2000 

Crescimento 

(%) 

Bom Fim 

11711 

11351 

-3.2 

Cidade Baixa 

19441 

16634 

-16.9 

Centro 

43252 

36862 

-17.3 

Santana 

23589 

21221 

-11.2 

Menino Deus 

30309 

29577 

-2.5 

Floresta 

25408 

14941 

-70.1 

Hípica 

7868 

10363 

24.1 

Guarujá 

2123 

2589 

18.0 

Serraria 

4528 

5775 

21.6 

Chapéu do 

Sol 

nc 

nc 

nc 

Ponta Grossa 

1711 

3290 

48.0 

Lajeado 

3942 

3425 

-15.1 


4. Conclusões 

A estimativa de áreas impermeáveis para 
o planejamento urbano é item fundamental aos 
gestores públicos pela correlação com a densidade 
habitacional conforme as referências citadas neste 
artigo. Relação esta diretamente ligada à infraestru- 
tura urbana como o saneamento, saúde e transporte 
já que retrata áreas de maior concentração popula- 
cional bem como tendências de avanço demográfico 
quando se comparam dados de censos realizados 
ao longo do tempo. 

Quanto à classificação realizada neste tra- 
balho, obteve-se inicialmente um erro da ordem 
de 40% na classificação da imagem de controle, o 
que conduziu pela análise da matriz de erro a um 
refinamento com maior número de amostras e a 
retirada de classes como água e sombras presentes 
na imagem. Deste modo, obteve-se um acerto global 
de 70% e um índice kappa de 0,57, considerando-se 
satisfatória a classificação. 

O estudo apresentado vem corroborar a 
tendência verificada e comprovada por estudos 
anteriores, demonstrando a aplicabilidade como 


ferramenta de planejamento urbano a estimativa 
de áreas impermeáveis como dado de entrada para 
modelos hidrológicos, bem como do ordenamento 
territorial quanto ao avanço de áreas urbanas e 
consequentes impactos ambientais associados à 
ocupação. 

Como conclusão deste estudo, houve um 
decrescimento populacional em regiões centrais, 
não havendo um impedimento para a densificação 
destas áreas pela população existente em direção 
às regiões periféricas, fato este, observado nacio- 
nalmente. 
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